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ABSTRACT 


The  scattering  of  E-polarized  cylindrical  electro¬ 
magnetic  waves  by  an  infinitely  long  dielectric  cylinder 
is  investigated.  For  small  incident  wavelength  the  slowly 
converging  series  solution  is  converted  to  a  sum  of 
integrals.  An  expansion  is  performed  so  that  the  integrals 
may  be  divided  into  two  classes  depending  on  the  relative 
locations  of  source  and  observation  points.  Those  in  the 
first  class  have  points  of  stationary  phase  and  an  asymp¬ 
totic  evaluation  gives  contributions  which  are  identified 
as  those  of  geometric  optics.  Ihe  remaining  integrals  arc 
evaluated  as  residue  sums.  A  plot  of  normalized  back- 
scattering  cross  section  vs.  normalized  radius  for  n  =  .4 
is  given. 


1. 


Introduction: 

The  scattering  of  electromagnetic  waves  by  an  infinite  dielectric 
cylinder  of  large  radius  has  been  treated  by  Beckmann  and  Franz  D]  .00 
among  others.  They  employ  the  Watson  transformation  to  convert  the  slowly 
convergent  aeries  solution  into  a  sum  of  Integrals,  Borne  of  which  are  evaluated 
by  the  method  of  stationary  phase  and  yield  geometric  optics  terms.  The  re¬ 
mainder  are  transformed  to  residue  sums  and  interpreted  as  diffraction  effects. 

The  pole  locations  they  use  in  evaluating  these  residues  are  only  qualitatively 
determined . 

In  the  following  the  authors  employ  the  Poisson  sum  formula  as  suggested 
by  Wu  [3]  to  obtain  integrals  similar  to  thoso  found  by  Beckmann  and  Franz. 
However,  the  present  work  is  an  extension  of  theirs  in  three  respects:  (l)}  wore 
precise  pole  locations  are  used  in  the  residue  computation;  (2)  ,  the  refractive 
index  of  the  cylinder  is  taken  to  be  less  than  one;  and  (3)^  numerical  results 
are  included. 


Formulation: 

If  on  electromagnetic  wave  emanates  from  an  infinite  line  Bource  with 
its  electric  vector  parallel  to  the  axis  of  an  infinite  dielectric  cylinder,  the 
governing  equation  in  the  coordinate  system  illustrated  in  figure  1  is 


where 


[Vr*e  +  G  (?,?')  ~ 


kM  *  1 


k, 


(i) 


k,j  k*  real. 


2. 


The  function  0  (  P,  r'  ),  together  with  ito  derivative  9r  Gi?  t  j  is 

continuous  at  the  cylinder  uurface  r»)rj«=  ix.  and  satisfies  the  radiation 


condition 
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The  solution  of  (l)  in  r  >  a.  ,  found  by  separation  of  variables,  Id 
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The  line  source  may  bo  at  either  the  primed  or  unpriced  point,  provided  r'  >  r 
For  large  cylinders,  i.o.  V,y  »  I  ,  the  series  (2a)  converges  very 
slowly,  and  it  in  convenient  to  employ  the  Poisson  sum  formula  CO  to  obtain 
o* 
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When  Debye's  asymptotic  forms  for  the  Honkel  and  Bessel  functions 
substituted  in  (3)  it  is  evident  that  for  0  <  8  ^  IT  every  integral  with 
has  many  stationary  phase  points.  These  contributions  may  be  separated  by 
expanding  R  ( V  )  (See  Appendix  A): 
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liubatltutlnn  ('0  in  (3),  interchanging  ouramatlon  and  integration,  and  carrying 
out  varioun  algebraic  mnnipulationn,  ve  obtain: 
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Kach  of  the  Integrals  (5a)  -  (5?)  which  has  a  ntationary  phase  point 
yieldB  a  character In  tie  ray  of  geometric  optica.  Some  of  theoo  rayo  can  reach 
only  a  restricted  region  which  depends  upon  the  relative  dielectric  constant 
of  the  cylinder  and  the  relative  position  of  Dource  arid  observation  points. 
With  the  index  of  refraction  N  =  y /X  !  and  the  observation 

point  at  )r  <  r 1  ,  three  regions  can  be  distinguished  (see  figure  l). 


Region  1  -  Forward  Scattering 

In  this,  tbe  "shadow"  region,  it  is  convenient  to  rewrite  GCr^rO  in 

the  form: 
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wliore  (  )  and  (be  )  an  well  au  (^il),  (‘pf),  and  (5b)  have  been  combined. 

The  denominator  of  (6c)  contains  the  factor  W^(y)H^fy)  which 

increases  exponentially  when  0  >  y  .  The  principle  contribution  of  those 
integrals  therefore  occurs  in  the  ranr.e  [~y  )•  .'Substitution  of  the  Debye 

asymptotic  foimn  in  thin  interval  yieldo 
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where 
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where  the  inverse  conine  functions  are  bounded  by  0  and  TT*  .  If  is  a 

point  of  stationary  phase,  the  condition  "9©/3w  =  Q  leads  to 
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This  equation  does  not  have  solutions  for  all  values  of  m  and  p 

For  a  fixed  value  of  p^  for  example,  the  left  side  of  (8)  has  a  maximum 
in  (-A,  A)  which  always  occurs  at 


\?r  ■  -A.  As  A  approaches  y  this  maximum 


6. 


value  approaches  2TTp  —XJ  ,  where 

'U  -  tos  -hco-s'ij/ktr)  —  geos'1  N/  . 


On  the  other  hand,  at  d<y  ~  A  the  value  is  minimum  and  approaches  "L7  as 
A  approaches  y.  Since  the  x'inbt  hand  bMg  of  (8)  must  be  within  these  limits, 
the  possible  values  of  vr\  satisfy 

2irp  —  (l_7  +  0)  >  2.Trm  >  \J  -  0  . 

But  m  lo  an  integer,  and 
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so  that 
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it  1b  clear  that  for  fn  *»  O  ,  (8)  has  a  solution  only  when  IT  -  The 
equation  of  the  line  separating  regions  2  and  3  io  9  >=  XT  (see  figure  l). 

The  stationary  phase  evaluation  of  (7a)  yields 
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where  the  subscript  0“  denotes  evaluation  at 

©„.=•  kJL]  -  2  [ciL- ( 'Jr/bj']  4 j 
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7. 


The  physical  Interpretation  of  these  results  an  the  rayu  of  geometric 
optica  io  piven  in  w  •  We  note  briefly  that  the  function  multiplying  kg. 
in  (%)  in  the  dintance  the  ray  travela  outside  the  cylinder,  while  the  multiplicand 

of  k,  corresponds  to  the  p  equal  paths  the  ray  traverses  inside.  Thus  the 

np-l 

rays  undei'go  p-  1  internal  reflections  an  iu  borne  out  by  the  factor  > 

since  the  quantities  Rj^  R-2tr  >  (  1  +  R|(r),  ^  (  1  +  T^,-,-)  arc  the 

Fresnel  reflection  and  trannminaion  coefficienta .  furthermore  when  m  <  p/j? 
the  ray  encircles  the  origin  m  timen  in  a  counterclockwise  direction,  while  for 
m  ">  p/2  1  the  clockwise  encirclements  are  Riven  by  p  -  m.  The  localization 

principle  M  an  well  an  Snell's  law  are  verified  by  the  interpretation  of  >V/k 
as  the  distance  from  the  ray  to  the  origin,  where  k  «  for  the  external  ray 
and  k,  for  the  internal.  If  of;  and  CC^aro  respectively  the  angles  of  incidence 
and  refract  Lon  we  have 
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which  la  Snell's  law.  These  quantities  are  illustrated  in  figure  2, 

The  Integrals  in  (6a)  and  (6b)  do  not  have  stationary  points  in  region  1; 
they  are  evaluated  an  residue  sums.  Only  the  lower  half  plane  solutions  of 
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are  used  to  locate  the  polos  of  R  (  1?  )  since  with  W<0  the  integrals  (6a) 
require  that  the  contour  be  closed  below  the  axis.  The  pole  locations  t 

as  computed  in  to  arc  shown  in  figure  The  integrals  in  (6b)  are  similarly 


H. 


evaluated  at  the  polen  of  rt  (  -  ^  )  in  the  upper  half  plane,  i.e.  at  =  Qy 

Tbuu  (6a)  in  equal  to 
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Under  the  condi  tion;!  -  %  a  Q  ( and  —y  -  Q  (j)  ,  the  auyinptotic  form 
of  thin  rent  due  i;um  in 
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p(XjiJ)  and  are  the  Uchttbo  polynomial:!  W,  and  Ai  (  Yj  )  in  the  Airy 
function  of  the  flrnt  kind,  Similarly  (6h)  becomen 
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where  the  derivative  ia  f'iven  by  (lib). 
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The  i.crtni!  (lla)  and  (12)  represent  Uie  well-known  creeping  waveu  bl  • 
which  rndiate  tangentially  a:;  they  propagate  around  the  cylinder.  The  exponential 
decay  depends  on  the  product  of  the  angle  and  •  Thun  the  only 

important  contributlonn  are  from  the  firot  encirclement,  i.e.  rvi  »  o  and  -  1. 
Physically,  the  waves  have  l-adlaled  away  do  much  energy  in  their  first  encircle¬ 
ment  that  their  next  contribution  is  negligible.  The  localization  principle  can 
be  applied  to  thooe  waves  also  since  ^  ^  X  for  the  counterclockwise  waves 
which  are  associated  with  rays  striking  the  cylinder  tangentially  in  {?  >  O  J 
for  the  clockwise  waves,  -  v) 9  X  —  X  ,  and  the  tangential  rays  are  incident 
in  &  <  O  .  The  geometrical  model  in  illustrated  in  figure  4. 

For  0  =. 'TT'  ,  the  terms  in  (lla)  and  ( 1 2 )  are  identical  since  the  geometry 
is  symmetrical  about  0*0  ,  T  .  Further,  the  residue  sum  for  m  -  o  does 
net  converge  when 

Cos  '(%/ kzr-)  +  cotf  ‘ (x/ktr> ) 

which  defines  the  common  boundary  of  regions  1  and  2.  Here  the  observation  point 
moves  into  the  region  of  direct  illumination,  and  the  Integral  with  m  »  0  gives 
optical  rays. 


Region  2 

In  region  2  it  is  convenient  to  rewrite  the  integral  in  (6b)  with  m  ■  0  as, 
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The  first  of  these  integrals  (13a)  is  evaluated  by  the  stationary  phase 


technique .  We  obtain  in  part 
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which  is  Juot  the  first  tom  in  the  asymptotic  aeries  for  (t/4-)  N0  ikd*- > 
the  free  apace  Green's  function  or  equivalently  the  direct  ray.  In  addition  (l3a) 
yields  the  singly  reflected  ray 
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and  is  the  solution  of 


cos' (fa-/ ks')  +  eos'ttir/ktr)  -  2cos''  -  0. 


Here  we  have  neglected  the  phase  contributed  by  R(  (\?)  for  *  >J 
since  It  is  of  a  lower  order  than  that  given  by  the  Hankel  functions. 

The  integral  (l3b)  iB  evaluated  as  a  residue  sum  at  the  poles  of  the 
integrand  in  the  upper  half  plane.  These  are  the  solutions  of  ~  O  , 


ag  well  as  the  solutions  of 
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The  aolutiong  of  (l4b)  in  tbo  upper  half  plane  are  ,  while  those  of  (l*+a) 

are  0^  and  (see  figure  5).  Thus  it  appears  that  double  poleB  exist  at 

vL  ;  however,  an  examination  of  equations  ( l4a)  and  (l4b)  Bbowe  that  these 

solutions  are  not  exactly  coincident  and  in  fact  depend  on  2  ^  Wi^CyJ  in 

this  part  of  the  Q  plane.  Dy  writing  for  the  roots  of  (l^a)  and 

for  those  of  (l4b)  we  show  in  Appendix  B  that  the  total  residue 
contribution  at  ^  is  negliRible. 

Thus  (l3b)  is  equal  to 
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12. 


It  is  possible  to  interpret  formulae  (l5a)  and  (15b)  in  termn  of  ray 
optics,  but  not  as  clearly  as  in  tbe  creeping  wave  situation.  Since  3g  2 ->f  + 
where  >  0  ,  the  localization  principle  relates  these  terms  to  tbe  rayn 

which  strike  the  cylinder  at  the  angle  of  total  reflection  above  the  line  9  =•  0  . 
The  firnt  phase  term  in  (l5Q)  is  the  product  of  ^  and  a  distance  consistent  with 
the  above  interpretation  (see  figure  6).  The  second  phase  term 
-  [ioi~'(fs  /kxy')  4-  cos~,{b9/kt.r)  -  Zcos~’  (\fs/%)  -  0 J 
is  approximately  equal  to 

iv>sL&  *-  cosT' (NA.fr')  +-  cos"'  (Nfi.fr)  —  2cos~'Nj 

J 

which  is  jVj  multiplying  the  angle  St  shown  in  figure  6.  Thus  this  part  of  the 
phaoe  ropreoents  a  wave  traveling  counterclockwise  around  the  cylinder  with 
propagation  constant  .  Since  Ti'vtCiPj)  >0  >  the  exponential  in  (15a)  has 

a  negative  real  part  and  the  wavo  decays  in  proportion  to  the  angle j  however, 
the  f?e.(  1)5  )  In  negative  and  so  the  physical  interpretation  is  not  completely 
satisfactory. 

Regioii_3 

The  re-nultn  derived  above  for  (13a)  and  (13b)  are  applicable  to 
observation  points  in  region  3-  The  integrals  (5c),  however,  as  was  mentioned  in 
connection  with  equation  (8),  have  no  stationary  phase  points  in  region  3«  After 
summation  (5e)  ia  equal  to 

00  irt  01  ,'OB 

_  i  f  NjW'j  Wktr)  e'  , 

3  Jjiy)  L6>)f  *  _  y  J/fy;  1  * 

J^Cy)  1  (16) 

which  is  evaluated  by  residues. 

For  9  >  O  ,  the  contour  may  only  be  closed  by  an  infinite  semicircle 
in  the  upper  half  plane.  Tbe  four  residue  series  obtained  from  the  poles 
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and  Y>,  (see  figures  3  and  5)  converge  qo  alowly  that  one  could  aa  veil  employ 
tho  aerloo  solution  (2  cl  ).  It  is,  however,  poaBible  to  approximate  (l6)  by 
deforming  the  path  of  integration  as  shown  in  figure  J.  The  integral  over  the 
new  path  in  noglocted  compared  to  the  residues  at  the  captured  poles,  -  v?5  , 

(see  Appendix  C),  Wo  note  that  the  rosiduea  at  the  captured  poleB,  —  , 

are  negligible  as  shown  in  Appendix  B.  Thus  (l6)  is  approximately  equal  to 

!*<■-  4'). _ i  '4  [a-  oiK-j  3 _ 


Ifxgh  _  y  &\ 

t*  L/W/yl  J )  1  ^ 


yL{ - I - 1 

V  i 

•  exp|cki^r'1'-  i\>sfki)L  +'f*r*'-Ws/kot-  -  z  J  op-  CVs/k.*)"*-'^  + 

—  i^s  [cc>$'[  Q$/kx.'r')  +  dos  1  (Os/ktrJ  —  2  t-os  '(iJs/x)  +-  Q  J  ^ 


(IT) 

where  the  summation  includes  only  captured  poles  and  the  derivative  is  given  by 
(15b).  In  region  3  the  inequality 

cas~'(NaJr‘)  -/  COS  ' ( Nojr)  —  2.cos'{N)  _>  0 

(18) 

is  satisfied  and  successive  terms  in  (r?)  have  decreasing  magnitudes. 

Tho  form  of  (17)  in  tho  same  ao  (15®)  with  the  sign  of  B  reversed; 
indeed  for  0=0  the  integrals  (l6)  and  (13b)  are  identical.  The  physical 
interpretation  of  the  residues  at  \)$  and  —  is  similar  except  that  the 

former  are  associated  with  rays  incident  in  0  >  O 

Plane  Waves  and  Back-Scattering 

To  specialize  the  above  results  to  the  case  of  an  incident  plane  wave  of 
unit  amplitude,  one  simply  applies  the  operator 

hm  -  4-i  (TTk^Yz/1  exp  [o  (rr/4  -  k,r‘)J 

to  the  previous  results ,  We  choose  here  to  consider  the  case  of  back-scattering 
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since  the  other  case  of  chief  interest,  the  scattering  cross  section,  cannot  be 
completely  treated  by  the  methods  developed  above. 

The  normalized  back- scat tering  crons  section  is  defined  by 


(teT  E5c*=o>|\.|elt|1 

r-f.oa  I 

where  11  is  the  ncattered  field,  and 

r< 


Er.  = 

■  r  f/,, w * k~x- *-) p" f  cos--‘-  V 

L P£oS«i  ,  | 

p  ■>  z  r*f  b.  i  ,  "  —l  J 

'  NJ  c.osot»- 


'  prtiZ** rj  ^  , 

coi  /  -a  r  J  ^ 


(19a) 


*  y  /i:)/jr/  V  *  expE^jtUn-Qrr  v- 

-8(1)  o  -  uJs  Or-  zcosWs/x)  -£%.}■ 


(19c) 


with  coo  oc^  ■  (  i  —  Po-*/;**-  )•■ 


’  and  cosoty  .  (  /  _  )s  , 


The  calculation  of  0Tq  was  programmed  for  the  IBM  7070  of  the  Brown 

University  Computing  Center.  The  geometric  optics  terms  in  (19a)  required  the 
iterative  solution  of 

sina^  »  N  sinoCj. 


and 


<*r  =  -  L<*i  -IT  (f/z  ~W)~J  . 
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The  creeping  wave  terao  represented  by  (19b)  were  found  to  be  negligible  and  only 
(19c)  contributed  to  the  diffraction  effect.  Figure  8  1b  a  plot  of  vs  kza. 

for  N  -  .4, 

The  reoultn  developed  here  do  not  apply  to  observation  pointB  near  6-1 ? 
and  the  bouiularleu  separating  regions  1,  2  and  3  since  these  lineB  are  caustics. 
Thus  the  forward  scatter  with  an  incident  plane  wave  has  not  been  completely 
determined  and  the  ncattcrlng  croon  ocction  cannot  be  calculated. 


Ajppendlx  A 

Tbe  power  oerien  (4a) 

OO 

"R , C v> J  -+-  C i  +- R,>C i  +-  RZP"' 

converges  by  the  ratio  test  when 

I  R*  I  <  I. 

On  the  real  v*  axiu  thin  condition  becomes  I  or  Rz  ^  I 

We  write 

r>  „  ni  E  ♦ 

so  that 


r  r*-=  !g!I±2ilIgjLz  pg» ) 

|Dl’-+  N‘  |E|*  _  N  +  Dfr  ) 


We  have,  however, 

Ft>*  OE*  -  (E*  P*  -f-  DB)  =  -  (0-P*XE -E*) 


where 


and 


D-  D* 


V)  / 

Hv>  (x) 


HlJ}  lx)  4-( 

m  wTiTywi3- 


* 


4-i 

> 


(A-l) 


(A-2) 


£  -  e 


vy  iH^iy)! * 


l6. 


so  that  (A-2)  is  positive  definite.  This  implies  that  (A-l)  is  less  than  unity 
and  the  series  converges. 


Appendix  B 

At  the  solutions  of  (lUa),  ,  the  residue  calculation  gives  for  the 

denominator  of  (l3h) 

H."V  J  x  H _  J/V(y)  1  1  l  hS’w  „  Hpiy)  l 

1  hvm  ^ 


3v?  H^cx; 


2t 

7T 


3  f  _h“w  _  v H,w(y?  7 
3^  1  TvJ7^;  Hi*;Cy>  J  4'  -> 


where  the  Wronsklan  relation  has  heon  employed.  Similarly  at  the  solutions  of  (lVb), 
V$j//  ,  we  have 

T  ,  y  ?  3  f  „/£%  wJ-^yn 

h_,  tp  j-v  cp  i  *wiryHpip> 3»i  -  *  .7^)  J 


2i  9  ] 

7r  9^1 


-  y 


(j) 


i. 


Adding  the  two  residues  and  setting 


IT 

2l 


3d  i  ~h£'V) 


1  { 

3d  l 


J-2CyJ  ) 

j-d  (y;  J 


J-sCy)  J  da" 


Vjl'  =*  4"  =  ^ 


the  result  is 


31,1  hSV  J  3»  i  Hj'w 


J-vi  - 


S\s 


4 


'X,  _ 


9. 

3d 


f_e_L  ? 
1  W’cyJ]1! 


j 3  r  Hpc*) 

^L*WT*> 


-  y 


H'fcp 

H^CyJ 


if 


2X*  c^~f  £46^  expCzritj) 

[hu’w]' 


17- 


which  is  very  small  since  Jr*f(Ojt)>0  and  ^  (j)  increases  ag  an 
exponential  function  of  ^ 

The  residues  of  (l6)  at  -  yield  the  same  resul^as  substitution 

of  —  0  for  \?  will  show. 

Appendix  C 

On  those  parts  of  the  path  shown  in  figure  7,  which  coincide  with  the  real 
axis,  the  integrand  of  (16) 

hSWo  hJW)  eil,e 


July)  vJ/l y)> 

Wc*J  J  Jo  iy)  J 

has  an  exponential  variation  given  by 


(C-l) 


Hj \ktf)  HyiM  Ll)e 


01, 


(C-2) 


When  liPl  >  x  and  ,  (C  -2)  is  very  Boall  and  decreases 

rapidly  with  increasing  |\?| 

The  exponential  variation  of  (c-l)  on  the  circular  part  of  the  path 
changes  from  (C-2)  on  the  arc  PQ  to 


e‘ 


OB- 


(C-3) 


on  the  arc  QRS  (see  figure  7)-  Since  (c*2)  is  symjietric  about  the  imaginary  axis 
and  [Hit}(y)J  1  is  very  Bmall  in  the  third  quadrant,  the  maximum  value 
of  (C-l)  occui's  on  the  arc  PQR,  Along  PQR,  (C-2)  increases  from 


exp  | [+ 2($Z- —  2.0  cosh'1  (0/x)  J 
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at  P,  to  exp  (  ii?CV  )  at  R.  However,  ( C - 3 )  muBt  be  employed  at  points  between 
Q  and  R,  and  the  additional  quantity  DtffyJ  varies  exponentially  from 
unity  at  Q  to  exp  (  -  i\>7F  )  at  R.  Thus  the  integrand  (C-l)  attains  its 

maximum  value  in  the  vicinity  of  Q. 

At  Q  the  exponential  variation  of  (C-l)  ia 

exp  2  [je*-- -v)  [to*1  (^^4. 


+  cos-1 


ztos"1  L'J/x) 


j 


and  the  equation  of  the  curve  passing  through  Q  is 


(C-t) 


Xw»Qy“-v*>,/*  -i)to*"t»/y)J  =0. 


(C-5) 

Since  a  general  examination  of  (C-4)  subject  to(C-5)  is  quite  complicated,  a 
numerical  example  will  be  considered.  Let  |<,r  ^  3,5"X.  N -yfx  «=  .5  y 

and  W|*s  Ii4x  .  An  approximate  solution  of  (C-5)yields  y}  =  /,+x  cxp^-iTr/4) 
and  (c-U)  is  approximately  equal  to  exp  !•  z7  ^e)]  From  ( X3 )  1  vliicii 

defines  region  3,  we  obtain  6  <  ,7£  radians  and  the  integral  over  the  semi¬ 
circle  is  negligible  compared  to  the  residues. 
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FIG. 6  PHYSICAL  INTERPRETATION  OF  RESIDUES,  EQUATION(l5a) 
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